Excitonic Aharonov-Bohm effect in a two-dimensional quantum ring 
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We study theoretically the optical properties of an exciton in a two-dimensional ring threaded 
by a magnetic flux. We model the quantum ring by a confining potential that can be continuously 
tuned from strictly one-dimensional to truly two-dimensional with finite radius-to-width ratio. We 
present an analytic solution of the problem when the electron-hole interaction is short-ranged. The 
oscillatory dependence of the oscillator strength as a function of the magnetic flux is attributed to 
the Aharonov-Bohm effect. The amplitude of the oscillations changes upon increasing the width 
of the quantum ring. We find that the Aharonov-Bohm oscillations of the ground state of the 
exciton decrease with increasing the width, but remarkably the amplitude remains finite down to 
radius-to-width ratios less than unity. We attribute this resilience of the excitonic oscillations to 
the non-siniply connectedness of our chosen confinement potential with its centrifugal core at the 
origin. 
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I. INTRODUCTION 

Recent advances in nanofabrication of quantum rings 
and dots by self-assemblingji"— lithographio^ or etching 
techniques^ have opened an active area of research both 
theoretical and experimental. In such systems electrons 
and holes arc confined in a small region and consequently 
the Coulomb interaction is enhanced. The existence of 
bound states of electron-hole pairs offers a unique oppor- 
tunity to explore the Aharonov-Bohm (AB) effect 
for excitons in quantum rings J^i^^ Despite the exciton 
being a neutral entity it has been predicted to be sen- 
sitive to a magnetic flux due to its finite size inside a 
quantum ringi ^^d^ In experiments, this sensitivity would 
show as an oscillatory dependence of both the optical 
transition energy as well as the oscillator strength upon 
the magnetic fiuxj^iii^ii^"— Theoretically, the excitonic 
AB effect has been studied by a variety of different ap- 
proaches. A short-range interaction between the electron 
and the hole has been investigated for one-dimensional 
(ID) ringS j^^d^d^d^ where also the effect of an external 
electric field can be included.— Intermediate models as- 
sume two-dimensional (2D) rings with narrow width un- 
der harmonic confinement and Coulomb-like interaction 
potentials between the electron and the hole^^ii^ or ra- 
dial polarized excitons when electrons and holes move in 
different circles.—"— The excitonic AB effect in 2D rings 
has been studied in models with harmonio^^ and geomet- 
j.| q15^27,28 confining potential or using a 2D attractive an- 
nular Hubbard model.— In all cases, the excitonic AB 
effect for neutral excitons has been argued to be sup- 
pressed in 2D as the width of the ring is incrcasedj^ Re- 
cently, experimental results in molecular-beam epitaxy 
grown nanorings made by AsBr^^ etching and on self- 
assembled InAs/GaAs quantum dots^ report oscillations 
in the binding energy of neutral excitons which may be 
accounted for by the excitonic AB effect. 

In this paper, we consider the excitonic AB effect in a 



confining potential that can be continuously tuned from 
strictly ID to truly 2D with finite radius-to-width ratio 
while preserving the central structure of a ring, namely, 
its non-simply connectedness due to an infinitely strong 
repulsion at the origin.— We present a simple analytic ap- 
proach to the excitonic problem when the electron-hole 
attraction is short-ranged . -'^'^'^^ We then study how the 
amplitude of the AB oscillations in the oscillator strength 
changes upon increasing the width of the ring. We find 
that the AB oscillations of the exciton ground state en- 
ergy decrease with increasing the width of the quantum 
ring, but nevertheless the effect remains noticeable down 
to regimes with radius-to- width ratios smaller than unity. 
This shows the robustness of the excitonic AB effect in 
2D. 



II. SINGLE PARTICLE STATES IN THE 
QUANTUM RING 

In the absence of Coulomb interaction, the Hamilto- 
nian of a single particle (electron or hole) subjected to a 
magnetic flux in a 2D quantum ring is given by 



= ^ (p. - q^Af + V{r, 



(1) 



where m^, pi and A are the effective mass, the momen- 
tum in the plane and magnetic vector potential, respec- 
tively. Here the subscript i = e,h refers to the electron 
and the hole, respectively. Electric charges are Qe = — e 
and qh = e. The quantum ring is modeled by an an- 
harmonic, axially symmetric potential with a centrifugal 
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The confining potential has the minimum at \ri\ = R (see 
Fig.[T]) and, for this reason, R will be used as a convenient 
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FIG. 1. Plot of the radial dependence of the confining poten- 
tial for 7 = 1 (short-dashed), 2 (dashed), 5 (long-dashed) and 
10 (solid). The strongly repulsive core at the origin remains 
very prominent even for decreasing 7. 



measure of the cfFective ring radius. Close to the mini- 
mum the potential reduces to the well-known displaced 
parabola, Vin) = 2Vb(r,/i? - 1)^ = {l/2)m,uj^{n- R)^, 
used in other theoretical studies of 2D quantum ringsi^ 
As — 7> 0, we see from ^ that the centrifugal core as- 
sures the survival of the essential feature of a ring: its 
repulsive barrier in the center. The effective width W 
of the quantum ring can be estimated from the single- 
particle ground state in the harmonic potential, namely 
W ^ (eo/2Vb)^/''i?, where eo = /2mR^ is the ring-size 
quantization energy)^ Notice that we assume that W is 
the same for electrons and holes, namely rrie = rrih = m. 
For the purpose of this work, all energies will be mea- 
sured in units of eo and we parametrize the strength 
of the confining potential by the radius-to-width ratio 
7 = R/W = (2Vb/eo)^^^- When 7 00, we approach 
the limit of a ID ring, whereas 7—^0 corresponds to an 
anti-dot geometryi^ Figure [1] shows the radial confining 
potential for different values of 7. 

In order to study the AB effect in the quantum ring, 
we choose A = (ArjAe) ~ (0, $ ft,/e27rr), corresponding 
to an infinitely thin magnetic flux piercing the plane of 
the ring perpendicularly. Here $ is the dimensionless flux 
through the ring and h/e the universal flux quantumi^ 
We note that due to the axial symmetry aroud the ring 
axis, all our results for energies have to be periodic in 
$ with period 1 and we hence restrict ourselves to the 
sector $ G [0, 1]. Then the Schrodinger equation for the 
electron in polar coordinates r^. =[re,0e) is written in 
dimensionless form 
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where Me = {neAe) represents the set of quantum num- 
bers for the electron, which are rig = 0,1,2,... and 



£e = 0,±l±2,... For brevity we define the dimensionless 
energy Ajv/e = Em^/^o and radial coordinate p^ = r^/R. 
The Schrodinger equation for the hole is the same aside 
from a change in sign in the linear term on $, and with 
a set of quantum numbers Mh = {rih, ih)- 

The normalized eigenfunctions of Eq. ([3]) are given by^ 
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where = V/e + 7^/4 with /e ^ — ^ defining an 
effective angular quantum number due to the confine- 
ment and the magnetic fiux. stands for the gener- 
alized Laguerre polynomials. The corresponding dimen- 
sionless energies are Xm^ = 7^(2ne + 1 + fee) — 7'*/2. The 
eigenfunctions and energies for the hole are the same as 
for the electron, with a effective angular quantum num- 
ber fh = £h + ^- The dimensionless zero point energy 
{ue = £e = ^ ^ 0) for the electron is Ag = 7^. 

Figure [2] shows the electron energy as a function of 
the parameter 7. We note that levels at higher values of 
the quantum number Ug become increasingly uncoupled 
for 7 > 5 and hence we expect to see nearly ID behav- 
ior for 7 values beyond this regime.— From the inset is 
clearly observed that in this 2D confinement regime, in 
the absence of interaction, the ground state energy for 
the electron (once sustracted the zero point energy) de- 
scribes an oscillation with the magnetic flux. For 7 > 5, 
the 2D oscillation is indistinguishable from the ID case. 
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FIG. 2. Dimensionless electron energy Xm^ as function of 
the radius-to-width ratio 7. Solid, dashed, dotted, dotted- 
dashed and double-dotted-dashed lines corresponds to Ue = 
0,1,2,3,4 and le = 0. For Tie = we also show le = 
5, 10, 15, 20. The flux is $ = in all cases. The inset shows 
the flux dependence of the energy of the ground state for dif- 
ferent 7. The black lines corresponds to ID results. 
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III. SOLUTION OF THE EXCITONIC CASE 

Within the effective-mass approximation, the Hamil- 
tonian of the interacting electron-hole pair is given by 
% = He + T^h + T^e-h, where T-Le-h is the interaction 
term. We model the excitonic interaction between the 
electron and the hole as a short-range potential of the 
form He-hire, rfi)/ Co = {2Tr)^/'^voRW6{re - Vh), where 
vq < parametrizes the attractive interaction strength. 
This contact interaction is the same used in Refs. [13 and 
[20I extended to a 2D case, where the area of the ring is 
2ttRW. In this definition we have carefully chosen the 
prefactors such that in the ID limit, 7 — > 00, the values 
of Vq become identical to the corresponding ID parame- 
ter and facilitate comparison with the results of Refs. [lil 
and[20l Thus, we express vq as — o/tt^ where a denotes 
the ratio of ID excitonic Bohr radius to ring circumfer- 
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Before continuing with the detailed study of the model, 
let us discuss some of the assumptions made and the lim- 
itations which we will encounter. Let us first emphasise 
that the restriction to equal electron and hole masses 
is simply a presentational convenience; all calculations 
shown here can easily be generalised to the case of un- 
equal massesSS but with a certain loss of clarity in the 
mathematical expressions. Nevertheless, we shall present 
some results for unequal masses later. The assumption 
of an infinitely thin current-carrying solenoid generating 
the magnetic flux $ is a theoretical construct. The ex- 
periments cited in the introduction all use a magnetic 
field B to generate the required $. This results in an ad- 
ditional, diamagnetic term proportional to B^, which we 
ignore here similarly to the experimental papers4ii'^'i^""— 
Certainly the most drastic assumption seems to be the 
(5-function potential for the two-particle interaction. Its 
use is of course motivated by our resulting ability to re- 
duce the computational difficulties as we will show below. 
Nevertheless, we wish to emphasise that there are also 
certain conceptual advantages associated with it: (i) in 
ID, the 5-function interacting many-particle problem has 
been solved exactly and hence the expression for the ex- 
citon binding energy on a line is known in terms of vq ^ 
(ii) In Ref. [ij, it was shown how the Bohr radius of the 
exciton similarly depends on vq. Both these parameters 
will of course vary when another form of interaction is 
considered. However, as also shown in Ref. [13, it is the 
ratio a introduced above which governs the strength of 
the AB oscillations. The effect of other two-particle inter- 
action potentials along the ring, when expressed in terms 
of a, will lead to similar AB oscillations and we expect 
at least qualitative agreement. Even for a long-range po- 
tential such as the Coulomb interaction, we expect this 
to hold as long as the overlap of wave packets on opposite 
sides of the ring, i.e. across the origin at r = 0, can be 
neglected. For the confining potential considered here, 
with its strong centrifugal core, this should be a rather 
good approximation. 

We construct the exciton eigenfunction as a linear com- 



bination of the electron and hole single-particle eigen- 
functions 

"iiVe^rh)^ ^ AA/^A/,^'0A/.(^e)'0Mh(T'/i) ■ (5) 

The Schrodinger equation for the electron-hole pair may 
now be cast in equivalent form 

Am^m^{Xm, + - A)'ipM,{,re)iiMhirh) 

M.Mh (6) 

+ {'i^f'^voRW5{re - rh)'i{re,rH) = , 

where A is the excitonic energy in units of eo- Fol- 
lowing an analogous procedure as in Ref. [10 the co- 
efficients Ahi^Mh are obtained multiplying Eq. ^ by 
"01/ ('''e)V'M, [rii) and integrating over the coordinates 
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G 



(7) 



Aj\/„ + Aj\/,, — A 
where we have defined 

Gm.m. = j d^rM'(r,r)^l,Jr)V'Ljr) . (8) 

Setting Te = rh = r in the expansion of (O, multiplying 
by V"!/' (^)V'm' (f) and integrating over the coordinates 
we finally obtain 

Gm'^M'i^ = ^ GM^MhPM^MhM'^M'^[A) , (9) 



with 



P, 



{2tt)^/'^VoRW 
AjWe + Am^ — A 



(10) 



To proceed we define the total angular momentum of 
the electron- hole pair in units of fi, as L = £e+ih - Because 
the system is axially symmetric, only states with L = L' 
can contribute to the excitonic system. This condition is 
even more restrictive under the dipole approximation, i.e. 
only excitons with total angular momentum L — Q can 
absorb light polarized perpendicular to the ring. There- 
fore Eq. ([To)) reduces to 
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X / drr7^M.(r)7^M,(r)7^A/,(r)7^M^(r) 
Jo 



(11) 



We note that in the limit 7 oo, the integrals in (fTTj) 
reduce to 1 / \/2ttRW for Ue = Uh = n'^ — n'f^ — and 
le + lh = ^e + ^h- Fo'^ other combinations of rie, n/i, n'^, and 
n^, the corresponding Pai^ai^m;.mi are less important 
due to the energy denominator in (1111) . 
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Similarly as for ID quantum ringSfi^i^ there is no an- 
alytical solution of Eq. © for finite values of Vq. In 
order to find approximate solutions, we hence need to 
cut off the sums at some maximally allowed values for 
Me and M/j. Figure [2] shows that for smaller values 
of 7, i.e. increasing ring width W, the level separation 
between the quantum states of the single particles is 
decreased. Therefore we use different €max and rimax 
values depending on our choice of 7. We have tested 
that our results do not change appreciably for the range 
of $ and Wo considered here. As in Ref. [13, Eq. © 
is reformulated as a standard left-eigenvalue equation 
Gk' = J2k GkPkk' (A) after mapping the quantum 
numbers Me,Mh K and M'^,M'f^ -)■ K' according to 

K ={(.+ ^max)(»^max + 1)^ + "-e("-max ^ I) + Uh + I SUCh 

that K, = 1, 2, . . . , (1 + 24iax)(l + J^max)'. The ex- 
citonic energies are obtained numerically by determining 
the values of A which result in the matrix Pkk' having 
an eigenvalue equal to 1. For a given A, all eigenstates 
can be found using ^ , ^ and ([5|) . An advantage of our 
approach is that it allows us to target the ground state 
directly by choosing a suitable starting value for A. 



IV. RESULTS 




FIG. 3. Exciton energy A at magnetic flux "1> = (solid lines) 
and 1/2 (dashed lines) plotted as a function of the interaction 
strength vq. For clarity, symbols are shown for $ = only 
and the results for 7 = 0.5 at $ = 1/2 have been suppressed. 
The thin dotted horizontal line denotes the onset of the single- 
particle continuum at $ = 0. The two thin black lines denote 
the ID limit for $ = 0, 1/2. The inset shows the amplitude of 
the excitonic AB oscillations as a function of the ratio m^jmh 
at interaction strength no ~ —2/-k'^. 



In Fig. [3] we plot the ground state energy A defined 
by Eq. (|9]) with ^^ax = 40 for 7 > 0.5 and rimax = 5 
(/fn,ax = 2916) for 7 < 5 or n„ax = 2 (/C^ax = 729) 
for 7 > 5 and as a function of uq for different values of 
7. For 7 = 0.5 we have used £i„ax = 30 and rimax = 6 
[Km&x = 2989). Here and in all following figures, when 
plotting the excitonic energies A, we have subtracted the 
zero-point energy 27^ of the non-interacting electron-hole 
system. 

We see that for all 7 and $ values, the increase of the 
interaction strength leads to the formation of a state 
with decreasing energy values below the onset of the free- 
particle continuum. Wc also compare in Fig. |3] the 2D 
exciton results with the ID ring studied in Refs. [13 and 
[20I . When the radius of the ring is 10 times its width 
(7 = 10) the 2D excitonic behavior is essentially indis- 
tinguishable from the ID results in the range of vq values 
studied. In particular, the differences between energies 
at different flux values at large 7 decrease. Nevertheless, 
for small 7^3 different magnetic flux values lead to 
quite distinct A values — even in a 2D quantum ring 
the exciton is sensitive to the magnetic flux. It is also 
interesting to note that for large 7, the bound state ener- 
gies are more negative for larger values of wq whereas 
for 7 < 1 we find evidence that smaller 7 values lead to 
smaller differences between different values of <&. 

The above results have been obtained assuming the 
same effective mass for the electron and the hole. To 
addrees the question of the robustness of the excitonic 
AB effect in a more realistic situation with different ef- 
fectives electron and hole masses, we have calculated the 
exciton energy as a function of the ratio me/m^. Let 



us define the amplitude of the excitonic AB oscillations 
as the difference of the exciton energy at $ = 1/2 and 
$ = 0, namely A(l/2) - A(0). The inset of Fig. Elshows 
this amplitude as a function of the ratio mg/m/i for dif- 
ferent values of the radius-to-width ratio 7 at interaction 
strength vq — — 2/7r^. In 2D rings (7 = 0.5) the en- 
ergy difference is almost constant and the assumption of 
equal masses is well justified. Upon approaching the ID 
limit, i.e. increasing 7, the amplitude of excitonic AB 
oscillations is reduced but the effect is still revealed. As 
an example, in common III-V compound semiconductors 
the ratio of the electron and ligh hole masses typically 
ranges from 0.6 to 0.9, and it can be seen in the inset of 
Fig. [3] that the reduction of the amplitude is small. 

Figure 3] shows the AB oscillations of the exciton en- 
ergy as a function of the magnetic flux $ within one flux 
period at different values of 7 and vq. In agreement with 
Fig. [21 we find that the AB oscillations are retained for 
radius-to-width ratios ranging from 7 = 0.5 to 10. This 
shows that the excitonic AB effect remains robust even 
in a ring of finite width. Upon increasing the 7 values 
for different $ values, we find mostly a moderate increase 
of the exciton energy, except in the vicinity of $ = 0.5 
where even the reverse tendency can be observed. 

In Fig. [S] we plot the amplitude of the excitonic 
AB oscillations for different interaction strength = 
— I/tt^, — 2/7r^, — S/tt^ as 7 is varied. We see that upon 
decreasing 7 from the nearly ID behaviour at 7 = 10 
towards 7 sa 1.5, there is only a slight decrease in the 
amplitude of the AB oscillations. Upon further decreas- 
ing 7, the oscillations weaken more rapidly, but even at 
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different $ and v g ■ This is similar to the ID behavior 
described in Refs.Il^ and[20l 

The oscillator strength, defined as 



F : 



I / d^r*(r,r) 



(12) 



is plotted in Fig. [T] Large values of F corresponds to a 
large transition matrix element from the exciton ground 
state into the vacuum. We find from Fig. [7] that the 
results for large 7 are in good agreement with the ID 
resultsJ^ And when decreasing the radius-to-width ratio 
7, the value of F does not suddenly drop to zero, again 
emphasizing the robustness of the excitonic AB effect in 
ring of finite width. 



FIG. 4. Exciton energy A as function of the magnetic flux 
$ for difi'erent values of interaction strength vo and radius- 
to-width ratio 7. The thin dotted horizontal line denotes the 
onset of the single-particle continuum at $ = 0. Only every 
second data point is shown for clarity in each curve. 

7 ~ 0.5, they retain about 30-40% of their original value. 
Results for other values of vq are similar. This again 
shows that even for rather wide rings, the excitonic AB 
oscillations persist in this 2D case. 
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FIG. 5. Amplitude of the AB oscillations A(l/2) - A(0) as 
a function of radius-to-width ratio 7 for different interaction 
strength vq. The dashed horizontal lines correspond to the 
ID limilji^, the vertical line denotes the 7 = 0.5 values. 

In Fig. |6] we plot the exciton probablity density 
|4'(re,r/i)p for different values of 7. We integrate j^'p 
over the radial coordinates pe, Ph and hence retain the an- 
gular dependence in Figs. I6alandl6c| whereas in Figs. [6b] 
and I6dl we integrate out the angular degrees of freedom 
and retain the pe, Ph dependence. From these figures we 
conclude that the exciton fills the available width of the 
ring. Figures |6a] and [6c] show that the exciton is indeed 
bound, i.e. the majority of the weight of j^'p resides along 
the diagonal 9^ = Oh- Analogous results are obtained for 



V. CONCLUSIONS 

Our results suggest that the excitonic AB effect orig- 
inally predicted for a ID mode l^^i^^ remains essentially 
unchanged when allowing for rings of finite widths as 
given by ((2|). We find that when we enlarge the ring 
width by one order of magnitude from I/7 = 1/10 to 
1, the magnitude of the AB oscillations drops by about 
15% only. In addition, we show that the qualitative be- 
havior of the oscillations both for the spectral position 
as well as the oscillator strengths of the exciton lumines- 
cence lines are again governed by the relative strength 
of attractive Coulomb interaction to ring radius. Our 
results are in good agreement with recent experimental 
observations where the magnitude of the excitonic AB 
oscillations was observed to be about 0.5 meV at binding 
energies of 4.35 mV for rings of about 11-22 nm radius 
and 7 « li^ We also note that our confining potential ^ 
has been chosen to retain its non-simply connectedness 
due to the infinitely repulsive centrifugal core at the cen- 
tre. Hence even for very wide rings, there is an essential 
difference with respect to the previously considered 2D 
confining potcntials i^^i^^i^^i^^" — This demonstrates that 
it is not so much the width or the exact shape of the 
confining potential, but rather the avoidance of the ring 
centre which is the important ingredient needed for the 
experimental observation of the excitonic AB effect. 

Last, we expect that the effects of external elec- 
tric fields^i^i^ and the formation of charged exci- 
tons remain similarly robust in 2D, whereas disorder ef- 
fecta21i22i2S should be less important than in the ID case. 
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FIG. 6. Dependence of the local wavefunction probability on (a,c) angular coordinates 6s, Oh and (b,d) radial coordinates 
pe = re/R, Ph = Th/R for 7 = 1 (a,b), and 10 (c,d) for wq = -2/7r^ and $ = 0. The values of |vl'(6le, 6*^)1^ and |*(pe,Ph)P have 
been normalized to lie in [0, 1]. The colors go from |*|^ e [0.9, 1] (red) to |*|^ G [0,0.1] (blue) in steps of 0.05 for (a,c), 0.2 
for (b) and 0.1 for (d). The thick black lines indicate in (b) and (d) the width W of the ring in each case. 




Work at Madrid was supported by MICINN (projects 
Mosaico and MAT2010-17180). 
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